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Abstract 

Adiabatic and configuration-fixed constraint triaxial relativistic mean field (RMF) approaches 
are developed for the first time and a new phenomenon, the existence of multi-chiral doublets 
(MxD), i.e., more than one pairs of chiral doublets bands in one single nucleus, is suggested for 
^'^^Rh based on the triaxial deformations together with their corresponding proton and neutron 
configurations. 
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The handedness or chirahty is a subject of general interests in molecular physics, ele- 
mentary particles, and optical physics. The occurrence of chirality in nuclear physics was 
suggested in 1997 Q and the predicted patterns ot spect.a exhibiting », .e.. the chital 
doublets bands, were experimentally observed in 2001 f2?]. 

Since the pioneer work on chirality in nuclear physics, lots of efforts have been made to 
understand the new Dhenomena and explore their possible existence in the nuclear chart, 

eg.. BQQBBaflQ . Experimentally, the chiral doublets bands have been identified in 

many nuclei in A ~ 130 mass region with the earlier suggested configuration 7ihii/2'^i'h~y^, 
A ~ 100 with TTgg^\ ® i^hii/2, and A ~ 190 with 7rhQ/2 ® i^'i'iy2- theoretical aspect, the 
chiral symmetry breaking was firstly predicted in the particle-rotor model (PRM) and tilted 
axis cranking (TAG) approach for triaxially deformed nuclei l| . It has been investigated later 
in hybrid Woods-Saxon and Nilsson model combined with shell correction method as well 
as the Skyrme-Hartree-Fock cranking approach [sj, and its selection rules for electromagnetic 
transitions have been discussed in a simple PRM j^. 

For triaxially deformed rotational nucleus, the collective angular momentum favors align- 
ment with the intermediate axis, which in this case has the largest moment of inertia. Mean- 
while, the valence particle and hole angular momentum vectors align along the nuclear short 
and long axis, respectively. These orientations maximize the overlap of the particle densities 
with the triaxial core and minimize the interaction energy. The three mutually perpendic- 
ular angular momenta can be arranged to form two systems with opposite chirality, a left- 
and a right-handedness. They are transformed into each other by the chiral operator which 
combines time reversal and spacial rotation of 180°, x = TTZ{7t). The breaking of the chiral 
symmetry in atomic nucleus is observed due to the quantum tunnelling between the systems 
with opposite chirality. 

The description for the quantum tunnelling of the chiral partners is beyond the mean field, 
as the usual cranking approach is a semiclassical model and the total angular momentum is 



not a good quantum number 



lll |. In contrast, the PRM is better suited for this Diirpose, 




which however is so far confined for one particle and one hole configurations only [Jj, Ill 
Its generalization for more particles and/or holes are still under development. Furthermore 
the deformation 7 and configurations in PRM are not self-consistent rather as inputs of the 
model. 

The relativistic mean field (RMF) theory has received wide attention due to its success in 
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describing the properties of nuclei and many nuclear phenomena for the past years |12 
It is interesting to search for nuclei with triaxial deformation and configurations with not 
only one particle and one hole but also multi-particle and multi-hole suitable for chirality in 
RMF theory. A multi-dimensional microscopic crankin g R MF model is very time consuming 
and has only been applied in magnetic rotation so far |l4i |. 

In this Letter, we will develop the adiabatic and configuration-fixed constraint triaxial 
RMF approaches to investigate the triaxial shape coexistence and the possible chiral doublets 
bands in A ~ 100 mass region. The existence of multi-chiral doublets (MxD), i.e., more than 
one pairs of chiral doublets bands in one single nucleus, will be suggested via the examining 
the deformation and the corresponding configurations. 

In RMF theory, the nuclei are characterized by an attactive scalar field S{r) and a 
repulsive vector field V{r) in the Dirac equation. 



{_,0! . V + V(r) + P[M + 5(r)]}^, = e,^, 



(1) 



which can be solved by expanding separately the upper and lower components of the spinor 
ipi in terms of eigenfunctions of the three-dimensional deformed oscillator in Cartesian co- 
ordinates 0Q,(r) and its time reversal state ^^(r) = T(f)a{r) with T = i&yK , 



a 



(2) 



In order to avoid the complex matrix diagonaliztion problem in Ref. [isj ]- 
as a product of three Hermite polynomials [3| , 



r) is written 



7! 



\ 



(3) 



The meson fields are expanded similarly with the same deformation /?o and 70 for the basis, 
but an oscillator length smaller by a factor of \/2 than that for the nucleons {bs = &o/\/2) 
has been used in order to simplify the calculations and to avoid the necessity of additional 
parameters. The oscillator frequency is hioo = 41^4"^/^ MeV. 

In order to check the convergence of the results with the number of expanded oscillator 
shells for fermions Ugf, the binding energies (upper) and the deformations (3 (middle) and 
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7 (lower) calculated with effective interaction PKl [l7] for ^'^^Rh as functions of nof are 
presented in Fig. [TJ The legends (filled circles, open circles, squares, triangle ups, triangle 
downs, and diamonds) represent the results obtained in basis with different deformation (/9o= 
0.0, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively). It shows that as long as Uof > 10, the binding 
energies as well as the deformations (3 and 7 obtained are almost the same. Furthermore 
they are independent on the deformation (3o of the basis. Similar convergence check for the 
bosons has also been done. In the following, a spherical basis with 12 major oscillator shells 
for fermions and 10 shells for bosons will be used, which gives an error less than 0.1% for 
the binding energy. 

In general the triaxial RMF calculation leads to only some local minima. In order to get 
the ground state for triaxial deformed nucleus, constraint calculations are necessary and in 
principal such calculations should be carried out in two-dimensional (3 and 7 plane. As such 
two-dimensional constraint calculations turn out to be expensive even for modern computer 
facilities, therefore the constraint calculations with (^Q\q + 2Q22j ■, i-^-, are carried out 
to search for the ground state for triaxially deformed nucleus. 

The energy surface and the deformation 7 as functions of deformation (3 in adiabatic con- 
straint triaxial RMF calculation with PKl for ^''^Rh are presented as open circles in Fig.|2fa) 
and (b) respectively. There are some irregularities in the energy surface. Furthermore some 
local minima are too obscure to be recognized and it is technically difficult to understand 
their corresponding single particle configurations. We therefore performed the configuration- 
fixed constraint calculation similar to what have been done in the non-relativistic case jisl ]. 
Starting from any point in the energy surface in the adiabatic constraint calculations, the 
configuration-fixed constraint calculation requires that the single particle orbits occupied 
are fixed during the constraint calculation, i.e., 

|(V^,-(/5 + 5/?)|V^,(/3))|^l. (4) 

The energy surfaces and the deformations 7 in configuration-fixed calculations with PKl for 
^°^Rh are given as solid lines in Fig. Efa) and (b) respectively. For each fixed configuration, 
the constraint calculation gives a continuous and smooth curve for the energy surface and 
the deformation 7 as a function of deformation (3. The irregularities in the adiabatic energy 
surface disappear. In comparison, the minima in the energy surfaces of the configuration- 
fixed constraint calculations become obvious, which are respectively represented by stars 
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and labelled as A, B, C, D, E, F and G. Their corresponding deformations (3 and 7 together 
with their binding energies are respectively given in Fig. Efa) and (b). It is interesting to 
note that, for each fixed configuration, the deformation 7 is approximately a constant (as in 
Figl2tb)), which means that the deformation 7 is mainly determined by its corresponding 
configuration. 

The energies for these minima including the ground state are within 1.3 MeV to each 
other but correspond to different deformations (3 and 7, which is a good example of the 
shape coexistence. The shape coexistence here is different from the spherical, oblate and 
prolate shape coexistence, e.g., in neutron-deficient Pt, Hg and Pb isotopes. It is the triaxial 
shape coexistence, i.e., for the ground state A: the binding energy E= 903.92 MeV (the data 
906.72 MeV), deformation (3 = 0.27 and 7 = 24.7 °, the excited minima B: E= 903.82 MeV, 
p = 0.25 and 7 = 23.3 °, C: E= 903.28 MeV, p = 0.30 and 7 = 22.9 °, and D: E= 902.69 
MeV, p = 0.22, 7 = 30.8 °. All the states A, B, C, and D have deformation /? and 7 
suitable for chirality . As these states are all in one single nucleus, if chiral doublets 

bands can be built on these states, it may lead to a new phenomenon, the existence of 
MxD in one single nucleus. Therefore, in the following, we will investigate their proton and 
neutron configurations in detail to see whether the particle and hole configurations required 
by chirality are available. 

Performing the configuration-fixed constraint calculations for the ground state, the single 
particle levels can be obtained as functions of deformation /3. The difference in single 
particle levels obtained by choosing other minima is negligible. The neutron and proton 
single particle levels obtained in such a way are presented in Fig. IHl The positive (negative) 
parity states are marked by solid (dashed) lines and the occupations corresponding to the 
minima in Fig. |21 are represented by filled circles (two particles) and stars (one particle). 
The corresponding quantum numbers for spherical case are labelled at the left side of the 
levels. As the energy curves become very stiff for small deformation in Fig. |2l the present 
configuration-fixed constraint calculations cannot be performed for deformation /3 < 0.06. 
Therefore we cannot distinguish the occupation of some low-j orbits, e.g., the last occupied 
neutron orbit for state A may come from 2^5/2 or 2^3/2, as marked in the figures. However, 
as the nuclear chirality is essentially determined by the high-j orbits, such illegibility does 
not influence the conclusions here. 

The binding energies, deformations P and 7 as well as the corresponding configurations 
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extracted from Fig. Elfor the minima A, B, C, and D in ^°^Rh obtained in the configuration- 
fixed constraint triaxial RMF calculations with PKl are listed in Table. HI Except the state 
D, in which there is no high-j neutron valence particle, we found the high-j proton and 
neutron configuration for the ground state (state A) as 7T{lgg/2)^^ ® ^^(1^11/2)^5 state B as 
7r(l(79/2)^^®z^(l/iii/2)"'^, and state C as vr(l5f9/2)~^®z^(l/iii/2)^, respectively. All of them have 
high-j proton holes and high-j neutron particles configurations, which together with their 
triaxial deformations favor the construction of the chiral doublets bands. It is interesting 
to note that the states A and B compete strongly with each other in energy. However, due 
to different parities, the states A and B do not mix up and are possible to produce the new 
phenomenon M^D. So far, a pair of chiral doublets bands with negative parity have been 
observed in ^°^Rh j^. It will be interesting to search for other chiral doublets bands in this 
nucleus. 

Similar to ^''^Rh, detailed constraint calculations are also performed for other isotopes 
in A ~ 100 mass region and the possibilities of M^D exist in other nuclei as well. The 
results will be published elsewhere in details. Here in Fig. IH the deformations (3 and 7 for 
ground states in ^^~^^^Rh, lo^-iie^^g ^^^^ loo-iisjj^ isotopes are presented. The shaded area 
represents the favorable deformation 7 for nuclear chirality. The nuclei, ^°^Rh and ^°^Rh, 
in which the chiral doublets bands have been observed P, 0| , are marked as filled circles. 
Apart from ^°'^Rh and ^°^Rh, the favorable deformation 7 for chirality has been found in 
102,108, iioj^]-^^ i08-ii2^g g^j^j ^^^In, which implicates that more chiral doublets bands can be 
expected in A ~ 100 mass region. For each isotope chain, the deformation 7 varies with the 
neutron number due to different occupation in i'hii/2- 

In summary, adiabatic and configuration-fixed constraint triaxial RMF approaches are 
developed for the first time to investigate the triaxial shape coexistence and possible chiral 
doublets bands. A new phenomenon, the existence of multi-chiral doublets (MxD), i.e., 
more than one pairs of chiral doublets bands in one single nucleus, is suggested in ^°^Rh 
by examining the deformation and the corresponding configurations. Similar investigations 
have also been done for Rh, Ag and In isotopes and also with other effective interactions. 
The observation of M^D is very promising in this mass region. 
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TABLE I: The binding energies, deformations /? and 7 and the corresponding eonfigurations for 
the minima A, B, C, and D in ^'^^Rh obtained in the configuration-fixed constraint triaxial RMF 
calculations with PKl. 



State 


E (MeV) 


(3 


7 


Configurations 


A 


903.9150 


0.270 


24.7° 


7r(l59/2)-'®i^{(l/iii/2)'[(24/2)' or (2^3/2)!]} 


B 


903.8196 


0.246 


23.3° 


7r(l59/2)~'<^i^{(l^ii/2)M(2(^5/2)' or (2^3/2)2] } 


C 


903.2790 


0.295 


22.9° 


7r(l59/2)"^ ® i^(l^ll/2)^ 


D 


902.6960 


0.215 


30.8° 


7r(l59/2)~^ ® ^^[(2^5/2)^ or (24/2)3] 
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FIG. 1: Total binding energy, deformation f3 and 7 calculated with PKl for ^''^Rh as functions 
the number of expanded oscillator shells for fermions ngf. The filled circles, open circles, squares, 
triangle ups, triangle downs, and diamonds represent the results as the deformation of basis /3o 
equals to 0.0, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively. 
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FIG. 2: (color online) The energy surfaces (a) and the deformations 7 (b) as functions of defor- 
mation P in adiabatic (open circles) and configuration-fixed (solid lines) constraint triaxial RMF 
calculation with PKl for ^^^Rh. The minima in the energy surfaces are represented as stars and 
labelled respectively as A, B, C, D, E, F and G. Their corresponding deformations /? and 7 together 
with their energies are respectively given in (a) and (b). 




FIG. 3: (color online) The neutron and proton single particle levels obtained in configuration-fixed 
constraint triaxial RMF calculations with PKl for ^''^Rh as functions of the deformation (3. Positive 
(negative) parity states are marked by solid (dashed) lines. The occupations corresponding to the 
minima in Fig. |21 are represented by filled circles (two particles) and stars (one particle) . 



9 



60 

40 

20 

t 
60 

40 

20 



O Rh isotopes 



Rh 



"Rh 



X 



O Ag isotopes 



Ag 



r 



60 

O In isotopes 

40 
20 



oAg 

1 02 . 1 04 . 



""■"'in 



112, 



p '" pep O 



0.1 



0.2 



0.3 



FIG. 4: (color online) The deformations /3 and 7 for the ground states in Rh (upper), Ag (middle) 
and In (lower) isotopes in constraint triaxial RMF calculations with PKl. The shaded area repre- 
sents the favorable deformation 7 for nuclear chirality. The nuclei ^''^Rh and ^'^^Rh, in which the 
chiral doublets bands have been observed, are marked as filled circles. 
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